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ABSTRACT 

We investigate if consequential angular momentum losses (CAML) or an intrinsic 
deformation of the donor star in CVs could increase the CV bounce period from the 
canonical theoretical value ^ 65 min to the observed value P„im ~ 77 min, and if a 
variation of these effects in a CV population could wash out the theoretically predicted 
accumulation of systems near the minimum period (the period spike). We are able to 
construct suitably mixed CV model populations that a statisticial test cannot rule out 
as the parent population of the observed CV sample. However, the goodness of fit is 
never convincing, and always slightly worse than for a simple, flat period distribution. 
Generally, the goodness of fit is much improved if all CVs are assumed to form at 
long orbital periods. The weighting suggested by King, Schenker & Hameury (2002) 
does not constitute an improvment if a realistically shaped input period distribution 
is used. 

Key words: binaries: close - stars: evolution - stars: mass-loss - novae, cataclysmic 
variables. 



1 INTRODUCTION 

Cataclysmic variables (CVs) are short-period binaries con- 
taining a white dwarf (WD) primary (with mass Mi) and 
a low mass main sequence secondary (with mass M2). The 
secondary fills its Roche lobe and transfers mass to the WD 
through the inner Lagrangian (Li) point. 

The main features of the orbital period distribution of 
CVs with hydrogen rich donors are the lack of systems in 
the 2-3 hr period range (the so-called period gap) and the 
sharp cut off of the distribution at around 77 minutes, as 
can be seen in Figure 1 (upper frame; e.g. Ritter & Kolb 
1998). 

So far theoretical models have been unable to repro- 
duce the precise position of the observed short-period cut- 
off and observed shape of the CV orbital period distribution 
near this cut-off. This is summarised in Figure 1. Systems 
that evolve under the influence of gravitational radiation 
(GR; Kraft et al. 1962) as the only sink of orbital angular 
momentum (AM) reach a minimum period at Pmin ^ 65 
minutes (Figurel, middle frame; Paczyriski 1971; Kolb & 
Baraffe 1999). The probability of finding a system within a 
given period range is proportional to the time taken to evolve 
through this region. We thus have 



NiP) oc ^, 
P 



(1) 



for the number N{P) of systems found within a given or- 
bital period range around P, and P is the secular period 
derivative at this period. We thus expect an accumulation 
of systems (a spike) at Pmin where P — (Figure 1, lower 



frame), while no such spike is present in the observed distri- 
bution (Figurel, upper frame). 

The orbital period evolution reflects the radius evolu- 
tion of the mass donor, which in turn is governed by two 
competing effects. Mass transfer perturbs thermal equilib- 
rium and expands the star. Thermal relaxation reestablishes 
thermal equilibrium and contracts the star back to its equi- 
librium radius. The minimum period occurs where the two 
corresponding time scales, the mass transfer time tM and 
the thermal (Kelvin-Helmholtz) time Ikh are about equal 
(e.g. Paczyhski 1971; King 1988). If tjv/ S> tKH then the star 
is able to contract in response to mass loss, but if Im <C Ikh 
the star will not shrink rapidly enough and will become over- 
sized for its mass. The position of the minimum period is 
therefore affected by the assumed mass transfer rate, and in 
particular by the assumed rate of orbital angular momentum 
(AM) losses. In this paper we investigate ways to increase 
the period minimum by increasing the mass transfer rate, 
and investigate ways to "hide" the spike by introducing a 
spread of Pmin values in the CV population. In particu- 
lar, we study the effect of a form of consequential AM loss 
(CAML) where the AM is lost as a consequence of the mass 
transferred from the secondary, i.e. Jcaml oc M2 (see e.g. 
Webbink 1985). 

In section 2.1 we outline our general model assumptions 
and introduce the prescription for CAML. In section 2.2 we 
present detailed calculations of the long-term evolution of 
CVs, and in section 3 we compare the observed short pe- 
riod CV period distribution with various theoretically syn- 



2 J. Barker, U. Kolb 



= 



10 



B 



E 

S- 6 



2 - 





-10 



>- -11 



-13 
0.05 



« 0.04 - 

E 

«! 

I 0.03 
o 

i 0.02 

o 
a 

^ 0.01 



T — I — 1 — I — I — r 



I ' ' ' ' I ' ' ' ' I 



I ' ' ' ' I 




I I I I I I I I M I I I I I I I M I I I I 
70 80 90 100 llO 




80 90 100 

Period (min) 

Figure 1. Upper frame: The observed period distribution of CVs 
with periods less than 116 minutes. Middle frame: Calculated evo- 
lutionary track in the orbital period versus mass transfer rate (M) 
plane. Lower frame: Period distribution expected from evolution- 
ary track in middle frame. 



thesized model distributions based on the calculations in 
section 2. 



2 THEORETICAL VERSUS OBSERVED 
MINIMUM PERIOD 

In this section we investigate possible solutions to the mis- 
match between the theoretical and observed minimum or- 
bital period in CVs. 



2.1 CAML description 

The orbital AM loss rate J of a CV can be written as the 
sum of two terms, 



J = Jsys + JcAML, 



(2) 



whore Jsys denotes the "systemic" AM loss rate, such 
as gravitational wave radiation, that is independent of mass 
transfer, while Jcaml is an explicit function of the mass 
transfer rate. We have 



dM2 
and 

Jcaml — » 



M2 -> 



(3) 



(4) 



Wc consider the general case in which the CAML mech- 
anism, along with nova mass ejections, causes a fraction of 
the transferred mass to leave the system. This fraction may 
be greater than unity as the primary may lose more mass 
during a nova outburst than was accreted since the last out- 
burst. 

We employ a generic prescription of the effect of a 
CAML mechanism, thus avoiding the need to specify its 
physical nature. Possible CAML mechanisms include a mag- 
netic propeller, i.e. a system containing a rapidly spinning 
magnetic WD where some of the transferred material gains 
angular momentum from the WD spin by interaction with 
the WD's magnetic field (see e.g. Wynn, King & Home 
1997), and an accretion disc wind (see e.g. Livio & Pringle 
1994). 

Our CAML prescription largely follows the notation of 
King & Kolb (1995). The AM is assumed to be lost via mass 
loss that is axis-symmetrical with respect to an axis A fixed 
at the WD centre but perpendicular to the orbital plane. 

We define a as the total fraction of mass lost from the 
secondary that leaves the system. We assume further that 
a fraction /3 (0 < /? < a) of the transferred mass leaves the 
system with some fraction / of the angular momentum it 
had on leaving the Li point. 

We also consider mass that is lost from the system via 
nova mass ejections, which over the long term can be consid- 
ered as an isotropic wind from the primary (see e.g. Kolb et 
al. 2001). This material will carry away the specific orbital 
angular momentum of the primary and will account for the 
fraction (a — /3) of the mass loss. We thus obtain 

Jcaml = ??& i^Mi + , (5) 

where we define rj — Pf as the CAML efficiency. For 
comparison with King & Kolb (1995) we equate this to 



M2 

Jcaml = v-j—J, > 0, 

M2 

and obtain 



(6) 



(7) 



For our calculations shown below we use the approxi- 
mation 

This is an adaptation of the expression given in Kopal 
(1959) and is accurate to within 1% over the range Q < q < 
0.4. 



2.2 Results of numerical experiments 

In this subsection we present calculations of the long-term 
evolution of CVs as they approach and evolve beyond the 
period minimum. For the computations we used the stellar 
code by Mazzitelh (1989), adapted to CVs by Kolb & Ritter 
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Figure 2. The increase in mass transfer rate and corresponding 
increase in minimum period for increasing CAML efficiency (from 
to 0.95 as indicated). 



Figure 3. Evolution of the stabihty factor D with mass ratio q 
for an initially marginally stable system, (Mi = 0.7 Mq, M2 = 
0.2Mq, t] = 1.0). 



(1992). Some of these evolutionary sequences are the basis 
for the theoretical CV period distributions we present in 
section 3 below. 



2.2.1 Consequential angular momentum loss 

We calculated the evolution of individual systems that are 
subject to CAML according to equations 6 and 7. We chose 
Ml = O.6M0 and initial donor mass M2 = 0.2Mq, with 
a range of CAML efficiencies < r? < 0.95 as shown in 
Figure 2. 

The systems initially evolve from longer periods towards 
the period bounce (right to left) at almost constant mass 
transfer rate. The minimum period increases with increasing 
CAML efficiency to a maximum of around 70 min for r) = 
0.95. 

Mass transfer stability sets an upper limit on the CAML 
efficiency. An obvious upper limit is 1, where all the angu- 
lar momentum of the transferred material is ejected from 
the system. Although the ejected material may carry more 
angular momentum than was transferred (as in the case of 
a propeller system where additional angular momentum is 
taken from the spin of the WD) this does not affect the net 
loss of orbital angular momentum. 

The maximum CAML efficiency still compatible with 
mass transfer stability could be smaller than unity. The sta- 
bility parameter D which enters the expression for steady- 
state mass transfer, equation 9 (e.g. King & Kolb 1995) must 
be greater than zero; this defines an upper limit on 77. 

-M2 = M2 (J^^ (9) 

A plot of D against q for an initially marginally stable 
system (Mi = O.7M0, M2{init) = O.2M0 and r] = 1.0) is 
given in Figure 3. The system initially exhibits cycles of high 
mass transfer rate M2 > 10~^MQ/yr {D close to 0) and very 
low meiss transfer rate M2 — > 0. The high states are short 
lived, on the order of 2 x 10^ years (see Figure 4) . The system 
finally stabilizes with D « 0.65. At around Pmin {q ~ 0.15) 
D starts to decrease further but always remains positive, 
settling at a value around 0.3. 




100 200 300 

age (10® yr) 



Figure 4. Mass transfer rate cycles for the initially marginally 
stable system as in Figure 3 

2.2.2 Structure of the secondary 

The tidal deformation of the secondary may have an 
effect on the period minimum. Calculations by Renvoize, 
Baraffe, Kolb & Rittcr (2002), [see also Kolb 2002] us- 
ing 3-dimensional SPH models suggest that the secondary 
is deformed in the non-spherical Roche lobe such that its 
volume-equivalent radius is around 1.06 times that of the 
same star in isolation. 

We mimic this effect in our 1-dimensional stellar struc- 
ture code by multiplying the calculated radius by a defor- 
mation factor A before the mass transfer rate is determined 
from the difference between the radius and the Roche lobe 
radius via 

-M2 = Mo exp (^-^^±^^^ . (10) 

Here Mo — lO~*M0yr~^ is the mass transfer rate of 
a binary in which the secondary just fills its Roche poten- 
tial and Hp is the photospheric pressure scale height of the 
secondary (sec e.g. Rittcr 1988). 

Figure 5 shows the effect on the minimum period and 
mass transfer rate for systems with various deformation fac- 
tors A, ranging from 1 (no deformation) to 1.24. The mass 
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Figure 5. Evolutionary track of a system with Mi = 0.6Mq and 
M2{init) = O.2M0 for various deformation factors (as indicated). 



transfer rate is seen to decrease with increasing deformation. 
This can be understood from the functional dependence on 
orbital period and donor mass in the usual quadrupole for- 
mula for the AM loss rate due to gravitational radiation 
(see e.g. Landau & Lifschitz 1958). Although the quadrupole 
formula is strictly valid only if both components are point 
masses, Rozzolla, Uryri & Yoshida (2001) found that the 
GR rate obtained using a full 3-dimensional representation 
of the donor star differs from the point-mass approximation 
by less than a few percent. 

It can be seen from the figure that with the deforma- 
tion factor 1.06 the minimum period increases from around 
65 min to around 69 min, consistent with Renvoize et al 
(2002) for geometrical effects alone. A deformation factor of 
around 1.18 was required to raise the minimum period to 
the observed value of ~ 77 min. This is somewhat larger 
than the intuitive expectation 



increase in radius : 



new period \ 
old period 



1.12(11) 



from Kepler's law and Roche geometry. 

In our calculations we consider the simple case in which 
only the geometrical deformation effects are taken into ac- 
count. The inclusion of the thermal effects considered by 
Renvoize et al (2002) have the likely effect of reducing Pmin, 
possibly by around 2% compared to the case with purely ge- 
ometrical effects 

One possible physical mechanism that could cause a de- 
formation factor above the value of 1.06 is magnetic pressure 
inside the star, as suggested by D'Antona (2000). 

We note that Patterson (2000) claims to find obser- 
vational evidence for "bloated" secondaries in short period 
CVs. On the basis of donor mass estimates from the ob- 
served superhump excess period he finds that the donors 
have 15 — 30% larger radii than predicted from 1 dimen- 
sional., non-deformed stellar models if gravitational radi- 
ation is the only AM sink. Even if true, this observation 
cannot distinguish between an intrinsic deformation of the 
donor star or the non-equilibrium caused by orbital AM 
losses in excess of the GR rate. 



3 PARENT DISTRIBUTIONS VERSUS 
OBSERVATIONS 

To test the statistical significance of the theoretically pre- 
dicted accumulation of systems near the period minimum 
("period spike") we calculated the period distributions of 
model populations for various assumptions about evolution- 
ary parameters. For each parameter a series of evolutionary 
tracks were generated, typically around 20. 

As systems evolve after the minimum period a point 
is reached (typically when M falls below lQ~^^{MQyr^^)) 
where numerical fluctuations in M become so large that the 
Henyey scheme no longer converges. The stellar code uses 
tables to interpolate/extrapolate the opacities and equation 
of state for each iteration, and in this region the extrapola^ 
tions become very uncertain. To extend the tracks we used 
a semi-analytical method as follows. 

The tracks were terminated at a value of log M2 = 
logM2(Pturn) — 0.3, where M2(Pturn) is the mass trans- 
fer rate at the minimum period for the track. The radius of 
the star for the final part of the track is approximated by 



R2 



RoM^, 



(12) 



where i?o and C, are assumed to be constant. The values 
of Rq and C, were determined from the final few data points 
for each track. (C takes a typical value of around 0.15 for sys- 
tems beyond the period bounce.) To generate the extension 
to the track we then calculated P from the Roche lobe con- 
dition, and M2 by assuming stationarity as in section 2.2.1 
(see Figure 1, middle frame for an example of an extended 
track). 

We weight the chances of observation to the brighter 
systems by assuming 



7 > 1-0. 



(13) 



for the detection probability. We tested the calculated 

model parent distributions for various values of the free pa- 
rameter 7 against the observed CV period distribution. A K- 
S (Kolmogorov-Smirnov) test is insensitive to the differences 
between the parent distributions. The greatest difference 
in the cumulative distribution functions (CDFs) of the ob- 
served and modelled distributions occur at the boundaries 
of the CDFs, i.e. in the least sensitive region for the K-S 
test (Press et al 1992). We thus decided to use the following 
modified test. 



3.1 The modified test 

For each parent distribution 10000 model samples each con- 
taining 134 systems were generated. ^ Each sample was 
tested against the model parent distribution using a x^ test, 
with 1, 2 and 4 minute bins. This range bridges the need for 
good resolution and significance of the x^ test which requires 
a minimum number of GVs per bin. The observed period dis- 
tribution was tested against the model parent distribution 

^ 134 was the number of observed CVs in the range 77 < 
P{min) < 116: Rittcr & Kolb 1998, internal update June 2001, 
as of July 2002 the number of systems in this period range is now 
152 though this does alter the values given by the test, the 
trends and hence the results remain unaltered 
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also, giving the reduced value Xobs- The fraction / of gen- 
erated samples with a reduced value less than Xo6s 
used as a measure of the significance level of rejecting the 
hypothesis that the observed distribution is drawn from the 
parent distribution. In the following we quote the rejection 
probability Pr=/. 



3.2 Magnetic and non-magnetic CVs 

Kolb & Baraffo (1999) noted that the observed distribu- 
tion of non-magnetic CVs (Figure 6, middle frame), and 
the observed distribution of magnetic CVs (Figure 6, lower 
frame) show no significant difference below the period gap. 
To test and quantify this we compared these distributions 
for P < 116 min, giving a reduced probability of 0.1213. 
Hence we cannot rule out that the distributions are drawn 
from the same underlying parent distribution. This is borne 
out by the results of comparing both distributions with a 
parent distribution that is flat in P (see also Table 1, en- 
tries F and G) which give similar rejection probabilities 
(Pr=0.709 and Pr=0.781, respectively). We thus find no 
significant difference between the two distributions. In the 
following we therefore test models against the combined 
magnetic and non-magnetic distribution of observed sys- 
tems. 

The lack of any distinct features in the combined ob- 
served period distribution (Figure 6, upper frame) does in- 
deed suggest an essentially flat distribution for the under- 
lying parent distribution. The flat distribution gives Pr = 
0.552 (for the 1 minute bin width, see Table 1). We use this 
value as a benchmark for the models discussed below. 



3.3 Parent populations 

We define a standard set of assumptions for simple parent 
population models as follows: 

(1) The primary mass in all systems is O.GMq. This is the 
value around which the majority of WDs in CVs are ex- 
pected to form (see e.g. de Kool 1992). 

(2) All systems form as CVs at orbital periods greater than 
2 hours. This is consistent with the secondary stars in CVs 
being somewhat evolved (see Baraffe & Kolb 2000). 

(3) The flux of systems through the period gap is constant. 
That is, sufficient time has elapsed since the formation of the 
Galaxy for a steady state to have been reached, so that the 
number of systems arriving at the lower edge of the period 
gap is just balanced by the number of new systems forming 
at orbital periods greater than two hours. 

(4) The CAML efficiency is set to 

(5) The systemic AM loss rate is Jsys = SJgr, and the 
deformation factor is A = 1.06, so that Pmin equals the ob- 
served Pmin = 77mm, thus enabling us to test the statistical 
significance of the spike. 

(6) Brightness selection factor 7 = 1.0 

A model population subject to these standard assumptions 
can be rejected with the probability Pr > 1 — 10^'*. 

In the following discussion of various population models 
we just quote any differences of individual models from this 
standard set of assumptions. 



10 I — I — I — 1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




ao 90 100 110 

Period (min) 



Figure 6. Observed orbital period distribution 76 < P(min) < 
116. Upper frame: all CVs; Middle frame: non-magnetic CVs; 
Lower frame: magnetic CVs. 

3.3.1 Age limit hypothesis 

It has been suggested that the currently observed short- 
period cut-off is not the true minimum period but purely 
an age effect (e.g. King & Schenker 2002). This would arise 
if systems that we currently observe have not had sufficient 
time to evolve to the true period bounce, as illustrated 
in Figure 7. Here systems at the currently observed short- 
period cut-off of around 77 minutes will continue to evolve 
to shorter periods for around another 8 x 10* years before 
reaching the period bounce (if Jsya = Jgr, A = 1.0). We ob- 
tain Pr = 0.659 for the period distribution generated from 
the single evolutionary track corresponding to Figure 7, cut 
at 77 minutes, (for a 1 mirmte binning, table 1, model B), 
quite close to the value for a flat distribution. In this period 
region AI ~ const (see Figure 1, middle frame). As P scales 
roughly as M, the discovery probability is roughly constant 
if 7 = 1. 

The same flat distribution would be obtained if some 
mechanism would cause systems to 'die' (e.g. become too 
faint to be detected) before reaching the period bounce. 
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Table 1. tests on various parent distributions. 
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^obs 


rejection probability 




rejection probability 


A 


1.003 


0.552 


0.728 


0.207 


0.738 


0.324 


B 


1.077 


0.659 


0.840 


0.332 


0.744 


0.329 


CI 


1.253 


0.831 


1.252 


0.754 


1.889 


0.947 


C3 


1.065 


0.609 


0.887 


0.380 


1.148 


0.669 


Dl 


1.379 


0.916 


1.422 


0.896 


1.991 


0.961 


D3 


1.122 


0.691 


0.913 


0.441 


0.938 


0.525 


El 


1.211 


0.786 


1.258 


0.779 


1.503 


0.848 


E3 


1.052 


0.580 


0.942 


0.454 


0.807 


0.361 


F 


1.107 


0.709 


1.339 


0.854 


1.514 


0.876 


G 


1.155 


0.781 


1.089 


0.636 


1.268 


0.753 



KEY 



A: Flat distribution versus total observed. 
B: Age limit versus total observed. 
CI; CAML versus total observed, 7 = 1. 
C3: CAML versus total observed, 7 = 3. 
Dl: Ml versus total observed, 7 = 1. 



D3: Ml versus total observed, 7 = 3. 

El: CAML plus Ml versus total observed, 7 = 1. 

E3: CAML plus Ml versus total observed, 7 = 3. 

F: Flat distribution versus magnetic CVs only. 

G: Flat distribution versus non-magnetic CVs only. 
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Figure 7. The age limit hypothesis. A system with current orbital 
period of 77 minutes will continue to evolve to the true period 
bounce at around 66 minutes for the next 8 x 10® years. 

Meyer and Meyer- Hofmcistcr (1999a) speculate that AM 
Her stars become propeller systems before the period 
bounce, and so are no longer observed as CVs as their accre- 
tion luminosity would be very low. For non-magnetic disc- 
accreting CVs Meyer and Meyer-Hofmeister (1999b) specu- 
lated that as the secondaries become degenerate, the mag- 
netic activity of the secondary reduces rapidly to zero. The 
disc would then be fed by non-magnetic material, thus re- 
ducing the viscosity of the disc plasma and vastly increasing 
the recurrence time. 

3.3.2 CAML efficiency and primary mass spectrum 

We now relax assumption (4) and allow systems to oc- 
cur with equal probability with any value of the CAML 
efficiency. This produces the period distribution in Figure 8 
(upper frame) for 7 = 1. A spike is still present, though now 



Table 2. test on CAML parent distribution versus total ob- 
served systems for various 7 



7 


v2 
^obs 


rejection probability 


1 


1.253 


0.831 


2 


1.167 


0.751 


3 


1.065 


0.609 


4 


1.106 


0.661 


6 


1.152 


0.957 



broadened, and peaked at around 87 minutes. The proba- 
bility distribution function (PDF) then falls with increasing 
period. This parent distribution gives Pr=0.831 (1 minute 
binning, see Table 1, model CI), somewhat larger than that 
for the flat distribution. 

The result of varying 7 is summarised in Table 2. At 
7 = 3 the rejection probability roaches a minimum value Pr 
= 0.61. This corresponds to the parent distribution shown 
in Figure 8 (lower frame). The broadened peak at around 
87 minutes is almost the same as for 7 = 1, but at longer P 
the PDF increases again, i.e. there is a minimum at around 
95 min. This is caused by a corresponding minimum of M 
along the tracks of Figure 2. 

So far we have assumed that all CVs in the population 
have the same WD mass. Observations (e.g. Ritter and Kolb 
1998) and population synthesis (e.g. de Kool 1992) show that 
a spread of WD masses is likely. To investigate the effect this 
has on the shape of the PDF near Pmin we relax assumption 
1 and adopt the WD mass spectrum calculated by de Kool 
(1992). 

The corresponding full parent distributions for 7 = 1 
and 7 = 3 give Pr = 0.916 and Pr = 0.691, respectively, for a 
1 minute binning (see Table 1 for full results) . These values 



The minimum period problem in cataclysmic variables 7 



n— I— p-i— I— I— I— p-i— I— I— I— p-i— I— I— I— p-i— I— rq 




period (min) 

Figure 8. Period distribution for a population based on a CAML 
spectrum (0 < »? < 0.95). Upper frame: 7 = 1.0. Lower frame: 
7 = 3.0. 
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Figure 9. Period distribution for a population based on a pri- 
mary mass spectrum (de Kool 1992) and CAML efficiency spec- 
trum (0 < r; < 0.95). Upper frame: 7=1. Lower frame: 7 = 3. 



are slightly higher (worse) than for the CAML efficiency 
spectrum population. 

If wc combine the effect of the primary mass distribu- 
tion and the CAML efficiency spectrum, i.e. relax assump- 
tions 4 and 1, we obtain the parent distributions shown in 
Figure 9. The PDF for 7 = 1 (upper frame) gives Pr = 0.786 
and exhibits a broad peak with a maximum at around 85 
minutes, followed by a gradual decrease with increasing pe- 
riod. The PDF for 7 = 3 Figure 9 (lower frame) gives Pr = 
0.580 and also shows a similar, though somewhat sharper, 
broad peak as for the case with constant WD mass. The 
values of Pr are lower (better) than either of the previous 
models alone (see Table 1 for full results). With 7 = 3 we 
approach a value similar to that of a flat distribution (0.552) . 

3.3.3 Deformation factor spectrum 

Here we we abandon assumption (5) and assume instead 

that the secondary stars are subject to various deformation 
factors A (as described previously in section 2.2.2). A min- 
imum value of A = 1.18 is used to set Pmin equal to the 
observed Pmin = 77min. Any A between 1.18 and a maxi- 
mum value \max has equal weight. The rejection probability 
for different Xmax are given in Table 3. There is a minimum 
in rejection probability (Pr=0.887) at Xmax — 1.35. The par- 
ent distribution generated for this value of Xmax is shown in 
Figure 10 and exhibits a gradually increasing PDF with in- 
creasing period to a peaJc at around 90 minutes, and then a 
gradual decrease. 

3.3.4 Initial secondary mass spectrum 

We now replace assumption (2) that all systems form with 
orbital periods greater than 2 hours with the other extreme: 
all CVs form with orbital periods of less than 2 hours. 



Table 3. test for the model based on a deformation factor 
spectrum (versus total observed systems, for 7 = 1.0) 



iiiaxiiiiuiii 
deformation factor 


y2 

■^obs 


rejection probability 


1.31 


1.392 


0.934 


1.33 


1.347 


0.915 


1.35 


1.297 


0.887 


1.37 


1.382 


0.928 


1.39 


1.452 


0.955 
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Figure 10. Period distribution for a population based on a de- 
formation factor spectrum. 
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Figure 11. Period distribution for a population based on an 
initial secondary mass spectrum. 



Table 4. test for the model based on a secondary mass spec- 
trum (0.13Mq < M2 < 0.17Mq) (versus total observed systems, 
for 7 = 1) 



binning (min) 




rejection probability 


1 


3.0878 


> 1 - 10-" 


2 


4.3233 


> 1 - 10"* 



Specifically, we assume that all CVs form with donor masses 
in the range O.ISMq < M2 < O.ITM© (this sets Pmax = 116 
min) , and that any M2 is equally likely. From this we obtain 
a parent distribution as in Figure 11. The PDF exhibits a 
sharp spike at the minimum period (here 78 minutes) and 
then a gradual decline with increasing period. The corre- 
sponding test results are given in Table 4. We conclude 
that if we were to include any secondary mass spectrum in 
the previous models we would effectively weight the PDF 
with a ramp function, similar to the one seen on the right of 
Figure 11. This would only increase the rejection probability. 



3.3.5 A contrived weighting? 

King, Schenker & Hameury (2002) constructed a (nearly) 
flat period distribution by superimposing individual ideal- 
ized PDFs with different bounce periods Pb according to 
a suitably tailored weighting. For the double box-shaped 
idealised PDFs modelled on the PDF shown in our Fig- 
ure 1 (lower frame) the required weighting is n{Pb) = 
exp[— 0.124(P6 — Po)] (Po is the observed minimum period). 
This weighting function effectively mirrors the shape of the 
sharply peaked individual PDFs. King et al. (2002) found 
that the range 78 < Pb < 93 is sufficient to wash out the 
period spike. It is clear that this procedure involves a cer- 
tain degree of fine-tuning for n{Pb) if the shape of the input 
PDF is given. Such a fine-tuning must surprise as the two 
functions involved presumably represent two very different 
physical effects. 

We applied the weighting n(Pj,) quoted in King et al. 
(2002) to our non-idealized model PDFs that involve the 
CAML efficiency and the deformation factor as a means 



Figure 12. Period distributions based on a deformation factor 
spectrum with 1.18 < A < 1.42 and n{P) = cxp[-0.124(Pi, - 
Po)] (solid line), n(P) = cxp[-0.07(P6 - Po)] (dashed line). The 
observed distribution (dotted line) is shown for comparison. 



to vary Pb. The weighting produced a marginally worse fit 

(Pr — 0.841 versus Pr — 0.831; 1 minute binning) for the 
CAML PDFs compared to the parent population based on a 
flat CAML efficiency spectrum we discussed earlier. In part 
this is due to the fact that the upper limit on ri does not 
allow a big enough range of Ph. In the case of the deforma- 
tion factor PDFs the fit marginally improved (Pr — 0.837 
versus Pr = 0.887; 1 minute binning, 1.18 < A < 1.35). It 
is possible to optimise the fit by adding systems with de- 
formation factors up to 1.42, and by using the weighting 
n{Pb) = exp[-0.07(Pi, - Po)], but this still gives the fairly 
large value Pr — 0.829 (sec also Figure 12). However, such a 
parent population is inconsistent with the observed distribu- 
tion for longer periods. As can be seen from Figure 5 systems 
that are subject to larger deformation factors would evolve 
into the period gap, hence the gap would be populated in 
this model. 

For completeness we show in Figure 13 the result of 
the superposition suggested by King, Schenker & Hameury 
(2002) if realistic rather than idealised PDFs are used. This 
model assumes additional systemic AM losses (5 — 11 x Jgr; 
no CAML, no deformation factor, 7 = 1) as the control 
parameter for varying Pb, and the weighting as in King et al. 
The pronounced feature just above 2 hrs orbital period is the 
result of the adiabatic reaction of the donor stars at turn-on 
of mass transfer (see e.g. Ritter & Kolb 1992). Such a feature 
is absent in the observed distribution. If deformation effects 
are taken into account the additional AM losses required to 
wash out the Pmin spike would cover a similar range but at a 
smaller magnitude. The resulting period distribution would 
be similar to the one shown in Figure 13 



4 DISCUSSION 

We have investigated mechanisms that could increase the 
bounce period for CVs from the canonical theoretical value 
~ 65 min to the observed value Pmin ~ 77 min, and ways 
to wash out the theoretically predicted accumulation of sys- 
tems near the minimum period (the period spike). Unlike 
King, Schenker & Hameury 2002 we focussed on effects other 
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Figure 13. Period distributions based on a spread of systemic 
AM losses 5 - 11 X Jgr and n{P) = e2;p[-0.124(Pi, - Pq] 

than increased systemic angular momentum (AM) losses, i.e. 
we assume that gravitational radiation is the only systemic 
sink of orbital AM. 

Wc find that even a maximal efficient consequential AM 
loss (CAML) mechanism cannot increase the bounce period 
sufficiently. As the real CV population is likely to comprise 
systems with a range of CAML efficiencies we would in any 
case expect to have a distribution of systems down to ~ 65 
min, rather than the observed sharp cut-off. 

We considered donor stars that are "bloated" due to in- 
trinsic effects, such as the tidal deformation found in 3-dim. 
SPH simulations of Roche-lobe filling stars. An implausibly 
large deformation factor of around 1.18 is needed to obtain 
a bounce period of ~ 77 min. 

A possible alternative identification of Pmin as an age 
limit rather than a period bounce (King & Schenker 2002) 
would limit the donor mass in any CV in a CV popu- 
lation dominated by hydrogen-rich, unevolved systems to 
> O.IMq. Any system with donor mass much less than this 
would either have an orbital period less than 78 minutes 
or would have already evolved beyond the period minimum. 
There are indeed systems with suspected AI2 < 0.1 AIq; good 
candidates are WZ Sge (M2 ~ 0.058; Patterson et al 1998) 
and OY Car (M2 ~ O.O7M0; Pratt et al. 1999). 

It is also possible that systems 'die' or fade before reach- 
ing the period bounce, and hence become undetectable as 
CVs. The fact that the very different groups of non-magnetic 
and magnetic CVs show almost identical values of Pmin (see 
Figure 6) strongly suggests that the physical cause for the 
potential fading would have to be rooted in the donor stars 
or the evolution rather than the accretion physics or emis- 
sion properties of the systems. 

Even if the bounce period problem is ignored we find in 
all synthesized model populations (except for the age limit 
model) a pronounced remaining feature due to the accumu- 
lation of systems near the bounce. We employ a modified 
test to measure the "goodness" of fit against the observed 
sample. An F-test (Press et al 1992) was also applied to the 
majority of 7 = 1 models and the same general trends ob- 
served. None of our syntliesised model populations fits as 
well as the distribution which is simply flat in orbital pe- 
riod (rejection probability Pr ^ 55%). Only models where 
brighter systems carry a far greater weight than expected in 



a simple magnitude-limited sample (selection factor oc M'' 
with 7 ~ 3 rather than ~ 1) achieve similar values for Pr- 
However, most of our models with 7=1 canot be rejected 
unambiguously on the basis of this test. 

Models designed to "wash out" the period spike by in- 
troducing a large spread of the CAML efficiency do generally 
better than population models based on donor stars that are 
subject to a large spread of intrinsic deformation factors. For 
all models the rejection probability decreeises if the full WD 
mass spectrum is taken into account, as this introduces an 
additional spread in the bounce period. Model populations 
where all CVs form at long orbital periods (chiefly above the 
period gap) give a much better fit than models that include 
newborn CVs with small donor mass. Adding these systems 
to the population introduces a general increase of the or- 
bital period distribution towards short periods, thus mak- 
ing the period spike more pronounced. This suggests that 
most CVs must have formed at long periods and evolved 
through the period gap to become short-period CVs. This 
is consistent with independent evidence that CV secondary 
stars axe somewhat evolved (Baraffe & Kolb 2000; Schenker 
et al. 2002; Thorstensen et al 2002). 

Recently, King, Schenker & Hameury (2002) con- 
structed a flat orbital period distribution by superimposing 
idealised PDFs that describe subpopulations of CVs with a 
fixed initial donor mass and initial WD mass, but different 
bounce periods. This superposition required a strongly de- 
clining rmmber of systems with increasing bounce periods. 
We repeated this experiment with a realistic PDF, but failed 
to obtain a markedly improved fit. 

In conclusion, we find that the period minimum prob- 
lem and the period spike problem remain an open issue. It 
is possible to construct CV model populations where the 
period spike is washed out sufficiently so that it cannot be 
ruled out unambiguously on the basis of an objective statis- 
ticial test against the observed CV period distribution. 
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